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Abstract 

The problem of interacting electrons moving under the influence of a 
strong magnetic field in two dimensions on a finite disk is reconsidered. 
First, the results of exact diagonalizations for up to N = 9 electrons for 
Coulomb as well as for a short-range interaction are used in the search 
for a peculiar ground state corresponding to filling factor 1/3. Not for 
the Coulomb, but only for the short-range interaction, can the 1/3— 
state be safely identified amongst the spectra of various filling factors 
close to 1/3. Second, the propositions of the concept of quasiparticles, 
as used in the hierarchical theory, are examined in view of the exact 
results for the disk geometry. Whereas the theory for the quasiholes 
is in complete accordance with the spectra, for the quasielectrons, fi- 
nite size corrections make an analysis difficult. For the quasielectron 
energy, an extrapolation to N — » oo is given and compared with the 
corresponding extrapolations of three different proposals for trial wave 
functions. While the limiting value for the best trial wave function is 
very close to the limit of the exact results, the behavior of the finite 
size corrections of the exact energies and of the trial wave functions, 
respectively, is qualitatively rather different. 
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1 Introduction 



The peculiar transport properties of a two-dimensional electronic system in a 
strong magnetic field, as seen in the integral quantum Hall effect (IQHE) ||T| 
and the fractional quantum Hall effect (FQHE) ||, are just a part of the rich 
structure of the phase diagram in the electron density vs. temperature plane 
@, |J. Much work was devoted in the last years to a study of the Wigner 
crystal phase at low densities |§ and, recently, to the interesting properties at 
filling factor 1/2 || |^, 0]. Here, we wish to reconsider the FQHE reporting 
work done on the few-particle problem. 

In both quantum Hall effects, there is a gap above the ground state at the 
magic filling factors in the ordered system. If a gap remains in the presence of 
disorder, a plateau in the Hall conductivity a xy is to be expected in the exper- 
iments (around filling factors v = p/q in the FQHE ; p,q - relative primes, q - 
odd); the longitudinal conductivity a xx vanishes wherever a xy is in a plateau 
region. Despite the experimental similarity in the transport properties, the 
microscopic models considered for the IQHE and FQHE are quite different. 
The IQHE is understood as a localization-delocalization transition, as shown 
by non-interacting electrons moving in a strong magnetic field under the ad- 
ditional influence of a moderate random potential ||. There, the energy gap 
is due to the cyclotron energy. In contrast, in the FQHE, even the starting 
point of a theoretical analysis, namely interacting spinpolarized electrons in 
the lowest Landau level without any disorder, is a non-tractable system itself 
because of the strong correlation. Here, even the existence of an energy gap 
in the many-particle spectra at the magic filling factors has to be proven. 
An important step towards an understanding was taken by Laughlin P, ITO, |TT| 
who showed the special character of the ground state at filling factors 1/q (q - 
odd) by presenting a trial wave function. This success was partially based on 
conclusions drawn from few particle calculations. Laughlin's notion of quasi- 
particles in connection with the hierarchical theory [|l!| |13| opened a way to 
explain the occurence of filling factors with numerators different from one, 
e. g. 2/5 . However, the hierarchical theory is rather qualitative and allows 
all filling factors with an odd denominator; thus, a further theoretical analysis 
becomes necessary. Moreover, a microscopic derivation of an effective quasipar- 
ticle Hamiltonian which would start from the original electronic Hamiltonian 
is still missing. 



A new and different point of view emerged in the work of Jain H, [L5L |16 



He circumvented the rather artificial construction of many particle states at 
v = p/q via quasiparticles which is done in the hierarchical scheme. Instead, 
he suggested explicit electronic wave functions for all observed filling factors, 
v = n/(2np + 1) (n,p - integers) and particle-hole transformed states. Jain's 
approach starts from n completely filled Landau levels, then multiplies the 
wave function with a symmetric Jastrow factor of power 2p ("adds 2p flux 
quanta per electron" ) in order to simulate the influence of the interaction and, 
finally, projects to the lowest Landau level. Numerical calculations for a small 
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number of particles and also the Fermi-liquid like behavior seen in the 
experiments || |7| near v = 1/2 corraborate Jain's approach. Nevertheless, 
the reason for its success remains unclear if one considers the quite different 
energy scales of non-interacting and interacting electrons, respectively JIT 



Other approaches attack the problem with field theoretical methods |19|, |20 



2jJ . These attempts describe already known features of the theory such as 
the Laughlin ground state, the nature of the quasiparticles and the collective 
excitation spectrum p2 . 



Both, the trial wave function approach and the field theories consider only 
points of the phase diagram of interacting electrons in the lowest Landau level, 
because the analysis focuses on definite filling factors. It should be emphasized 
that the approaches considered so far have another thing in common: they do 
not depend very much on the specific form of the electron-electron interaction. 
On the one hand, this can be seen as expressing a universality of the FQHE. 
On the other hand, one would like to check such an assertion. A nice concept 
for a discussion of the interaction is that of the pseudopotential coefficients, i. 
e. the eigenenergies of the two-particle problem, introduced by Haldane p3 . 
In the lowest Landau level, the interaction is completely determined by these 
coefficients. The picture emerged, that the largest of these determines the 
FQHE ground state, not to be reached from the non-interacting ground state 
by perturbation theory, while the other coefficients can be included perturba- 
tively and do not lead to a big difference. Is this true? Below, we will show 
the numerical spectra for two different interactions (i) the full Coulomb inter- 
action and (ii) a pseudopotential interaction, where only the first coefficient is 
non-zero. 

At first glance, the choice of an appropriate geometry for studying such a sys- 
tem seems to be only a matter of convenience. Three geometries were studied 
in the past: the disk geometry with a symmetric gauge and open boundary 



conditions |L0|, the torus geometry with periodic boundary conditions based 



on the Landau gauge, cf. pi], and the spherical geometry where the electrons 



move in the constant magnetic field of a magnetic monopole in the centre [13 



The analytical theory of Laughlin |Tl[ was formulated on the disk and also nu- 



merical work was done for this geometry [T(J, |26], |27|. Soon, the advantages 
of the spherical geometry with an additional symmetry and without bound- 
ary attracted a lot of interest. It became the favorable choice allowing exact 
diagonalizations for up to ten electrons at filling factor 1/3 P8| and the study 
of the hierarchical theory, because the lower energy levels are well separated 
from the rest of the spectrum. 

In this paper, we come back to the disk geometry for several good reasons: 

- From the experimental point of view, one should study a planar geometry, 
particularly, if one wants to account for boundaries and contacts. 

- The topological differences between the different geometries can very well 
change the results. E. g. , the ground state at v = 1/q is nondegenerate in 
the disk and spherical geometry, whereby in the torus geometry it is g-fold 
degenerate 



3 



- The Laughlin theory was originally formulated for the plane. Then, in the 
current understanding of the FQHE, the plateaus of the Hall conductivity in 
the vicinity of the magic filling factors are a result of the disorder. Up to now, 
calculations including disorder study the energy gap between the ground state 
and the excited states, see e. g. [^2], |33|, [34] ; the calculation of the transport 
properties of a disordered and interacting electronic system in a magnetic field 
is an outstanding problem. Still, current theories of disorder use a planar ge- 
ometry, see [[35)1 . 

- If one wants to verify theories stressing the importance of edge excitations 
in the FQHE |36 |, one needs a geometry with an edge. 



In this work, we present the results of a numerical study of the low-lying 
eigenstates of the FQHE Hamiltonian at and in the vicinity of filling factor 
1/3 for the disk geometry. We compare the Coulomb with a pseudopotential 
interaction. While on the sphere, the differences between both interactions 
are not so marked, on the disk, the pseudopotential interaction plays a special 
role and this will be shown below. With our results, we evaluate properties 
of quasiparticles and thus check the foundations of the hierarchical theory, by 
comparing the exact spectra for small numbers of quasiparticles with energy 
expectation values of trial wave functions. 

The paper is organized in the following way: in Section 2, the model with the 
interactions is introduced. Section 3 describes shortly the numerical method 
used for the diagonalization. In Section 4, we set out to identify the stable 
state in a system with a finite number of electrons, and we compare the result 
with Laughlin's ground state wave function. The notion of quasiparticles is 
discussed in general in Section 5. Then, the results for the quasihole are dis- 
cussed in Section 6, those for the quasielectron in Section 7. Here, we have 
extended our recent analysis |37j to a larger number of electrons. Section 8 
gives the summary. 



2 Model 

We consider electrons of charge — e (e > 0) moving in the x-y-plane under 
the influence of a perpendicular, constant magnetic field B = —Be z (B > 0). 
They interact via an isotropic, translationally invariant interaction V(|r*— r'|). 
We study the disk geometry and then, the appropriate gauge is the symmetric 
gauge A = f (y, -x,0). 

The one-particle problem is equivalent to an isotropic two-dimensional har- 
monic oscillator, whose Hamiltonian Hq and angular momentum L can be 
expressed by two pairs of commuting Bose operators a,a^ and b,tf, so that 
H = hcu c (tfb + |) with uj c = 1^ and L = h(a)a — Wb). Then, the one-particle 
basis is 

\n : m >= [n!(n + m)!]-5 (fet)n ( a t)*+m | > ^ 

with n = 0, 1, . . . ; m = —n, . . . , 0, . . .. The energies E n = Tiuj c {n + |) are de- 
generate in each of the Landau levels n with respect to the angular momentum 
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m. Because the gap between adjacent Landau levels grows linearly with B, 
whereas the coupling constant for a 1/r interaction increases only with the 
square root of B, in the strong field limit (FQHE), the Landau level mixing 
can be neglected and the Hilbert space of the one-particle basis (0) is restricted 
to n = 0. Analogously, if the g-factor in the Zeeman term is not too small, the 
electrons are spin polarized (about the necessity to include the spin, see p8|). 
The basis (|l|) becomes in real-space representation for the lowest Landau level 

<p m (z) = [2n2 m m\]--i z m e~^\ (2) 



where z = x + iy and all lengths are given in units of l c = yh/\eB\. 
We want to calculate the energy spectrum for a finite number N of electrons 
at a given filling factor. The "filling factor" z/, naively is the ratio of the 
number of electrons to the number of available one-particle states. In the 
disk geometry with background potential, this latter number is given by the 
number N$ of flux quanta - through the area of the background. For the 
neutralizing background potential, we take a superposition of all one-particle 
states with an angular momentum from m = to m max ; this leads in the 
thermodynamic limit to a homogeneous background charge density. Then, 
outside the disk with the area 2tiN$ = 2ir(m max + 1), the background charge 
density drops quickly to zero. Thus, there are just m max + l one-particle states 
with a maximum expectation value of the area not exceeding the area of the 
background. Therefore, v' — -S- — — ^—rr- For the numerical diagonalization, 
we restrict the single electron Hilbert space to only h + 1 one-particle states 
with m = 0,1, ... ,h. The Hamiltonian to be diagonalized numerically then 
reads: 

1 h 

H — — 2-j Wmi,m2,m3,m4 C mi C m2 C m3 C mi 

mi ,m2 ,m3 ,m,4=0 

h TTLmax Z/'^ m,max 

m=0 m'=0 m,m'=0 

where the kinetic energy term has been already dropped since it is a constant 
for fixed N. The c m , are Fermi annihilation and creation operators and the 
two-particle matrix element is 

WW.ma.ms.n* = / / <P* mi ( Z )fm 2 W) V{ \z ~ z' \ )<f m3 (^> m4 (z) d?Z(Pz'. (4) 



In the following, we specify our two different choices for interaction potential 
V(|;z|) and confinement. The complete model (§) for the case of a Coulomb 
interaction with background, the Coulomb model (CM), is the first one. The 
second is a short-range interaction, the short-range model (SRM). In contrast 
to the first case, here is no neutralizing background and the confinement is 
due to the truncation of the one-particle basis at h — m max . In this case, we 
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N-l _ N-l 



will use a slightly different definition for the filling factor, v — , . — 
because of reasons which will become clear below. The difference between v 
and v' disappears in the thermodynamic limit. Then, in the case of the SRM, 
the Hamiltonian reduces to the first term in (H) with h = m max . 
In first quantization, the first term of H can be decomposed as follows 

TV oo 
i>j=l m=0 

is the projector which annihilates in a two-particle wave function all 
components except the one with relative angular momentum m. This decom- 
position in angular momenta was originally introduced by Haldane for the 
sphere. The V m are the eigenvalues of the two-particle problem and are called 
pseudopotential coefficients. Because of the fermionic nature of the problem, 
only odd m (m — 1, 3, . . .) contribute to the sum in (H). The W TOl|TO2 , m3 ,m4 are 
expressed by the V m as follows : 



4 _ m\+m,2 

" "„ X 

=1 m=0 



W m H(2 mi m i !)-5 £ {m 1 + m 2 )\m\V n , 



oo /__ \ / \ oo 




The Coulomb model is specified by the pseudopotential coefficients 

v cm = VK(2m)\ 

2 2m + 1 (m!) 2 ' 1 ' 

Here and below, all energies are measured in units of the coupling constant 

2 

of a Coulomb potential. The SRM is specified by the coefficients 



ySRM = \/^ V £™ = 0. (8) 

Here, only electrons with relative angular momentum one interact and repel 
each other. An interaction in real space, which yields Vj^ RM as pseudopoten- 
tial coefficients in the lowest Landau level, was introduced by Pokrovsky and 
Talapov f3~9fl : 



N 

H 1 st M = V 1 SRM £ V 2 5 2 (r!-r}), (9) 
i>j=i 

hence the name short range model. Why should one study a SRM ? The 
ground state of @ at the magic filling factors v, where the experiment shows 
the FQHE, cannot be reached by a simple perturbation theory from non- 
interacting electrons. If, on the other hand, the non-perturbative ground state 
of a SRM is already very close to the true ground state of ([J), then the higher 
coefficients (V m> i in our case) can be switched on perturbatively and do not 
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change the essence of the results. In this SRM should be considered as 

a generic model for the FQHE. Indeed, the Laughlin wave function is the exact 
ground state for v = 1/q for a SRM with V m > q = [|13|, |40| . Particularly 



on the sphere, this reasoning is very successful. In addition, Haldane showed 
for v = 1/3 that if the V m >3 become too large compared with V\, then the 
preference of the Laughlin wave function as the ground state wave function is 
destroyed [p3| . 

For future reference, we give the translation operator. The generator of the 
magnetic translations is t = p — eA and commutes with the Hamiltonian for 
an infinite system. For N electrons, the magnetic translation is a product of 
unitary one-particle operators 

N N i N 

f ? = J] if = exp(* £ &) = exp(-^= ]T(a+r - 0,0) , (10) 
i=l i=l v * j =1 

where we have used again the bosonic operators a, a + . 
The action of the operator on the basis (Q) yields 

Tf VmOzi) = [2n2 m m\]- l 2 {z x - 0™ e-^ 2+ ^''\^ 2 . (11) 

Next, the factors {z\— S,) m and e^ 1 ^* on the r. h. s. of (|TT|) can be expanded with 
respect to ^, £*. Thus, the transformed function ( a "ring" of radius \j2(m + 1) 
around is a superposition of all the degenerate eigenfunctions of the lowest 
Landau level with angular momenta reaching up to infinity. For an electronic 
system of restricted size, one needs to truncate the angular momentum at 



3 Numerical treatment 

We want to study the energy spectrum, particularly the ground state and 
the lowest excited states, of if at a fixed filling factor. Then, the electron 
number iV determines the area of the system expressed by the degeneracy 
N$> or by m max . The dimension of the fermionic many-particle Hilbert space 
(^tv) = ( mm ^ +1 ) grows dramatically with increasing N. However, also for a 
finite system, the total angular momentum M = J2m=o mc L c m is conserved. 
Thus, the diagonalization can be performed in a Hilbert sub-space with fixed 
M, i. e. , in the "block M". The allowed values of M are M min = N(N -l)/2 < 
M < M max = N(2m max -N+ l)/2. 

For the diagonalization of (§) in a block M, we construct the many- particle 
basis in terms of the Slater determinants of the one-particle states. Then, 
the sum of the (1 < i < N, < < m max ) is M and all one-particle 
angular momenta have to be different (fermionic case). The dimensions of the 
blocks M are symmetric between M m i n and M max . The matrix with M being 
in the middle has the largest dimension and its size limits the feasibility of the 
diagonalization. The number of these many-particle basis states with fixed M, 
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gp(Nip, N, M), can be determined either numerically by a recursive formula or 
from the following generating function [41 



x 



N 

,N(N-l)/2 TT 
1=1 



1 — X 



1 — X' 



(JV*-(1-1)) 



M=M, 



g F (N ,N,M)x 



(12) 



As an example, the largest Hilbert sub-space for v — 1/3 and 10 electrons has 
M = 135 and the dimension 246 448. The maximum dimension is the same as 
that in the spherical geometry for the same number of electrons at the same 
filling factor, cf. @. 

The W mi ,m2,m 3 ,m4 i n ©h determined by the pseudopotential coefficients via 
(^D, are used to calculate numerically the matrix elements between the many- 
particle states; the diagonal elements are sums of the difference between the 
direct and the exchange terms, the off-diagonal elements are only single dif- 
ferences. With increasing N, more and more of the matrix elements are zero 
because of the angular momentum conservation and, therefore, the matrix 
becomes sparse. At filling factor v — 1/3, for N = 6, the matrix M = 45 
(dimension 338) contains 6108 non-zero off-diagonal matrix elements out of 
56953, whereas for N = 9 and M = 108 (matrix dimension 45207), there are 
3016844 non-zero matrix elements of 1021813821, i. e. only 0.3% off-diagonal 
elements are non-zero. The number of non-zero off-diagonal elements grows 
approximately as the matrix dimension to the power 1.25. We have checked 
that there is no further (accidental) decomposition of a block M into subblocks 
as long as we use the above basis of Slater determinants. 
For the diagonalization, we apply two methods. For N < 6, we calculate all 
energy levels by combining the Householder and the QR method. For electron 
numbers N > 7 when the matrices become quite large and sparsely occupied, 
we use a Lanczos method in which we restrict the search for energy levels to 
the low-energy region. 

In this work, we concentrate on the filling factor 1/3 and calculate the ex- 
act energy levels. We do not study overlaps between eigenfunctions and trial 
wave functions because our conclusions can be drawn already on the basis 
of quantum numbers and expectation values. In the spherical geometry, in 
contrast, one needs also the eigenfunctions for a determination of the angular 
momentum L 2 . 

4 The stable state 

In the following, we describe the results of our numerical calculations for a 
filling factor of precisely 1/3, where the experiment shows the FQHE. We 
study first the CM, then the SRM. 

How can one identify in a series of spectra for various m max at fixed N the 
one at filling factor 1/ql An important hint is given by the trial wave function 
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for filling factor 1/q, the Laughlin wave function fT0|| : 



N 



. . . ,z N ) = T[ (zi - Zj) q evp 

q 

i>j=i 



1 N 

_ I ~ |2 

4 ^ 

^ i=l 



2; 



(13) 



Its total angular momentum is M = qN(N — l)/2 = M^(q) ; we will abbrevi- 
ate Mjv(3) in the following as M^. (|T^) exhibits a large overlap with the true, 
numerically determined ground states for small N in the CM. This was shown 
only for up to four electrons on the disk JL(| and for the corresponding wave 
function on the sphere for up to ten electrons |^S] . There are many arguments 



why (|13D is so successful. In particular, if q is not too large (q < 70) [Q, it 
describes an incompressible, homogeneous liquid 0, For a finite number 



of interacting electrons on the sphere, the requirement of a homogeneous and 
nondegenerate state makes the identification of special spectra corresponding 
to filling factor 1/q simple, because the ground state must have zero total an- 
gular momentum jP| . 



We want to identify - in the disk geometry - in a series of spectra for the two 
interactions, CM and SRM, the ground state of the filling factor at 1/3 as a 
special state which we call the stable state and which we compare with the 
Laughlin wave function. There will be problems in the case of the CM, but 
the identification will be straightforward for the SRM. 

We start with the CM (with background potential included) For an 

overview, the complete spectrum of a system with 6 electrons on the disk 
calculated by the exact diagonalization procedure is shown in Fig. la ordered 
with respect to M. We have chosen m max = 17 so that the naive filling factor 
is v' = 1/3. We use h = m max in (|3|); numerical tests show that larger h 
have very small influence on the energy values but just increase the dimension 
of the Hilbert space. The ground state is nondegenerate and has M 6 = 45 
as predicted by the Laughlin function. Still, it is not a priori clear that the 
ground state for this choice of m max is a finite N approximation of the 1/3- 
state. We will investigate this now by changing m max ; generally, as we increase 
m max from its minimum value N — 1 (at v' = 1), the spectrum moves to larger 
M and the ground state angular momentum increases. In order to show this 
in detail, we compare for 7 electrons the low-energy levels in the vicinity of 
v' = 1/3 where m max varies between 16 and 21 (Fig. 2a - f); the spectrum at 
p 1 — 1/3 is given in Fig. 2e. At first glance, the most promising candidate for 
a particularly stable ground state seems to be the one in Fig. 2a, m max = 16, 
where the gap between the ground state energy level and the first excited state 
is the largest. But this m max does not correspond to v' = 1/3. Further, the 
state with one flux quantum more (m max = 17) does not show any signature of 
a degeneracy in the ground state in contradiction to the quasiparticle picture 
which is discussed in detail below in Section 5. Finally, the ground state spec- 
trum with the formally correct v' = 1/3 given in Fig. 2e does not exhibit the 
expected total angular momentum for seven electrons Mj = 63, but instead, 
M = 57. Table 1 summarizes the dilemma. It gives, from 5 to 7 electrons, the 
total angular momentum M g of the exact ground state for a range of filling 
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factors around 1/3 expressed by N$ = m max + 1. This range is chosen such 
that within it, the ground state angular mommentum M g increases through its 
Laughlin value M^. Also, m max and the angular momentum of (|T3| ) are given. 
The table shows that for the CM and N = 5,6, 7, an identification of the 1/3- 
ground state cannot be brought into agreement with what the Laughlin theory 
predicts. Either the maximum single particle angular momentum m max or the 
total angular mommentum M g does not agree with flT3|). Certainly, finite size 
corrections could in principle, in the limit N — > oo, remedy this problem. But 
the tendency seems to be in the opposite direction. The situation does not 
improve at all for eight electrons where the ground state at v' = 1/3 has an 
angular momentum 77 unequal to Mg = 84. 

Disregarding this difficulty, we extrapolate the ground state energy per parti- 
cle at v 1 = 1/3 for the block with Laughlin's value M = from data for up 
to 9 electrons to the thermodynamic limit. In the case v' = 1, the finite size 
corrections scale as -7=. Assuming this behavior also at u' — 1/3 yields 

E (-) 1 1 

= -0.409510 + 0.036220^ + O(-). (14) 



Fig. 3 shows the exact data and the extrapolation. Eq. ( |14|) is in good agree- 
ment with various other exact diagonalizations and with energy expectation 



value calculations using Laughlin's trial wave function on the disk 27] and on 



the sphere p8|, P| p| 



We now change the model in two steps to arrive at the SRM and do a sim- 
ilar study there. At first, starting from the spectrum Fig. la, we omit the 
background. This changes the shape of the spectrum dramatically, see Fig. 
lb. Without background, the repelling electrons tend to occupy the largest 
available angular momenta. In such a case, the Laughlin wave function cannot 
describe the ground state correctly. In the second step, changing the interac- 
tion to the SRM, we calculate the spectrum without background still using the 
same m max = 17. Then, the ground state of this spectum is highly degenerate, 
see Fig. lc. This shows that a change from the CM with background to the 
SRM using the same external parameters (N,m max ) does not necessarily yield 
a nondegenerate ground state, as we expect at the stable state. 
In order to find the nondegenerate ground state, we decrease m max successively 
and find for N = 6 at m max = 15 really a nondegenerate state of zero energy 
(Fig. le, for comparison see also the Figures 4a-f for 7 electrons including the 
stable state in Fig. 4c). For completeness, we show in Fig. Id also the spectrum 
for the case of the SRM with a background where only Vi-parts contribute to 
both electron-background and background-background interaction in (|3|). 
The nondegenerate ground state as seen in the spectra of Fig. le and 4c for 
6 and 7 electrons, respectively, can be immediately identified as the Laughlin 
wave function ( |T3| ) for N = 6 and N = 7. It is well known that ( |TB| ) is the 
exact and only solution of zero energy in the SRM at v = 1/3 because it is 
the only iV-electron fermionic wave function not containing a component with 



relative angular momentum one for any two electrons [13, EM HQ). The spectra 
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display the expected ground state angular momentum M = (in Fig. le, 
45 = Mq and in Fig. 4c, 63 = M7). The situation is similar in the case of the 
sphere where the stable state has L 2 = and is thus nondegenerate and homo- 
geneous. Eventually, this finding justifies our use of v which differs from the 
naive definition v' by a finite size correction, since v — 1/3 for the parameters 
of Fig. le and 4c. 

After having compared the spectra of the CM and the SRM at filling factor 
1/3 we conclude that the SRM seems to be - on the disk - the most promising 
model for a discussion of the concept of a stable 1/g-state and of quasipar- 
ticles. This is true although thermodynamic instabilities are expected in this 
model [ fjEQ . In contrast to the case of the spherical geometry, where the results 
for the CM do not differ qualitatively from those for the SRM |23|, |47[, the 



SRM suggests itself for a study in the disk geometry. Thus, from now on, we 
will exclusively consider the SRM. 

As a side remark, we supply an argument supporting the nondegeneracy of 
(|T3|) in the disk geometry. We try to construct a state degenerate with flT3"|) by 
applying a magnetic translation (commuting with H for the infinite system) 
to ( |T3"D and show that this leads back to (|I~3|). 

N 

f^i( Zl , ...,z N ) = . . . , z N )e-^\ 2 J] e W (15) 

9 q . 

i=l 

An expansion of the product w.r.t. £* yields eigenfunctions with total angular 
momentum M^(q) + 1 

00 1 

EEEJWO pg(^...,^)*i(^---,^v). (i6) 

1=0 k=0 j q 

The pjpj are symmetric polynomials of degree /, k counts the maximum power 

of any variable zi (0 < k < I) occuring in , and j enumerates the individual 
polynomial for fixed I and k. If is an exact eigenstate, as it is e. g. for 

the SRM, then for an infinite large disk all the functions Pu ^i are again 

1 

degenerate eigenfunctions, but with v = q ^~^ +k ■ However, we search for 



eigenfunctions with v — 1/q. Thus, we have to put k = which leaves in ( |i~6"D 
only the term I = 0, the Laughlin state (|T3"D. This shows how a translation of 
(|13|) followed by a restriction to the orignal area of the system (filling factor) 
leaves (|i~3D invariant (see for the case of other topologies |3l|). A similar 
procedure will prove useful in the case of the quasiparticles see below. 
Finally, we comment on the degeneracy in the spectra of systems without 
background potential for small M. M-blocks between M min and M min — (N — 
1) + frtmax < Mjsr(q) are entirely unaffected by the confinement (this is only 
true without background). Energy levels with different M are degenerate as 
the result of the center of mass motion (in Fig. lb and lc the part of the spectra 
between M — 15 and M = 27, and in Fig. le the part between M — 15 and 
M = 25). In some sense, even the Laughlin wave function is also such a state 
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because enlargement of the system does not influence the energy of the state 
with zero energy. 

It should be mentioned that we have refrained from calculating the overlap 
between Laughlin's wave function and the numerically calculated ground states 
because the agreement in M, in connection with the energy of the state being 
zero, suffices for our conclusion. Summarizing our attempts to identify the 
1/g-Laughlin state in numerically calculated spectra we see that for the SRM, 
the ground state can be identified uniquely already for a finite system. In the 
CM case, this is not possible. While on the sphere, the spectra of the two 
interactions for finite N are supposed to be related by perturbation theory, on 
the disk, they look completely different. 



5 Quasiparticles 

After having identified the stable state at the filling factor 1/3 in the SRM, 
we turn to the ground state properties at filling factors nearby. Laughlin 



Til was the first who introduced the notion of quasiparticles at 1/q in the 
disk geometry. Here, we look for a definition applicable in our finite system. 
Starting from the stable state, we create the two kinds of quasiparticles by 
increasing and decreasing, respectively, the number of flux quanta through the 
area. The quasihole is then defined as the bulk ground state for a system with 
one additional flux quantum through the area covered by the electrons, i. e. 
m max = q(N — 1) + 1, and a quasielectron is the bulk ground state for a system 
with one flux quantum less, i. e. , m max = q(N — 1) — 1. This leads to a smaller 
and larger filling factor, respectively, with a deviation of the order 1/N from 
1/q. 

In general, the quasiparticles are fundamentally different from ordinary quasi- 
particles of a Fermi liquid in that they cannot be constructed from the original 
fermions by a process of switching on the interaction adiabatically. They have 
various descriptions as anyons [12], bosons jl3| or fermions ||1 1|| . They are not 



elementary particles as the electrons, but charge deviations from the homoge- 
neous density of the Laughlin state. There are two applications of this notion 
of quasiparticles. The first is the description of excited states at an unchanged 
filling factor 1/q. In this picture, quasihole and quasielectron, separated by 
an infinite distance, form a quasiexciton [Q which is a short -wavelength- 
excitation {k — > oo) [^, |30], fOJ. The energy of this excitation is the sum of 



the two quasiparticle energies; this is the gap seen in the activation measure- 
ments of p xx [Q. Otherwise, for small k, a collective theory yields the low 



lying excitations f22 



Here, we are interested in the other use of the quasiparticles in which they de- 
scribe the ground state properties near 1/q and in which they are the elements 
of the hierarchical theory flTBf . This theory sets out to understand the oc- 



curence of filling factors p/q (p ^ l,q- odd). Starting with N free electrons in 
N$> energetically degenerate states in the lowest Landau level and switching on 
an interaction (e. g. the SRM-interaction), one gets a stable state for the spe- 
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cial filling factors v = 1/q, i. e. for iV$ — 1 = q(N — 1). The hierarchical theory 
now draws an analogy from the electrons to the quasiparticles. These can oc- 
cupy of the order of N degenerate states. Creating a macroscopic number N' of 
quasiparticles one gets in the presence of a quasiparticle-quasiparticle interac- 
tion a new stable state, the "daughter state on the first level of the hierarchy" , 
for the special filling factor 1/p of the quasiparticles, i. e. for N = p(N' — 1) (p 

- even for bosonic quasiparticles). Then, the electronic filling factor is given 
by 

N-l N^oo 1 

V (N — l)q ± N' ~* q±Y 1 ] 

For q = 3 and p = 2 this gives filling factors of 2/5 and 2/7. Iterating the 
scheme, one can find successively all filling factors between and 1 with odd 
denominators in an unique manner on different levels of the hierarchy. 
We want to emphasize that there are several propositions which the quasipar- 
ticles have to fulfill for a justification of this theory. 

- The lowest one-quasiparticle levels must be energetically degenerate at least 
in the thermodynamic limit, and there must be an energy gap to the other 
states of higher energy. 

- For a two-quasiparticle system (two flux quanta more or less), it must be 
possible to identify the low energy states as two free quasiparticles plus an 
interaction contribution. I. e. , these states must be separated from the ener- 
getically higher states which means that the interaction must be small enough 
to leave a gap between the two-quasiparticle states and the neglected upper 
states. 

- The low-energy levels for more than two quasiparticles, particularly for a 
macroscopic number of them, must be describable by an effective Hamiltonian 
with up to two-particle interactions as determined from the two-quasiparticle 
spectra. Higher interactions are not allowed. This is the strongest condition 
and the test for the validity of the hierarchical theory itself. 

In summary, the low-energy behavior of the effective quasiparticle Hamiltonian 
has to reproduce the low-energy behavior of the original fermionic problem. 
For the simplest non-trivial case of quasiholes at v = 1, this has been success- 
fully checked on the disk with a special interaction |51| . 



In the past, there were a lot of efforts to justify quantitatively the quasiparticle 
picture. Almost all authors used the spherical geometry and thus neglected the 
effect of boundaries. The work was mainly concentrated on the determination 
of the quasiparticle energies at various filling factors [IT], [27J and then, conclu- 
sions were drawn for the hierarchical theory. Beran and Morf [38] extracted the 
effective interaction of quasielectrons on the sphere and estimated the ground 



state energy and the gap of the 2/5-stable state. Endesfelder and Terzidis |52 



extrapolated the interaction of quasiholes and compared the spectra of a small 
number of bosons interacting via this interaction with electronic spectra of the 
same number of quasiholes. There is a large asymmetry between quasiholes 
and quasielectrons p3| . Their energies are of a different order of magnitude, 
and the subspaces of low lying excitations have different dimensions for the 
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quasiholes and the quasielectrons. Thus, it was tried to explain these different 
dimensions with a difference in the quasiparticle-interaction: both quasiparti- 
cles are bosons, but quasielectrons have a hard-core interaction. ]54j . 
In the two following sections, we want to study whether a degeneracy of the 
quasiparticles is found in the disk geometry. In the spherical geometry, this 
degeneracy holds trivially because the angular momentum of the quasipar- 
ticles is 1? = + 1) and, therefore, there are N + 1 degenerate states 
(//-degeneracy) |j30| . In contrast to the authors cited above, we investigate 
the quasiparticles in the disk geometry, where the situation is more difficult as 
a result of the boundary. 



6 Quasiholes 

A quasihole is created by adding one flux quantum to the system at v — 1/q 
by increasing the area of the stable state system by one flux quantum, i. e. , 
m max = (N — l)q + 1. For the short-range model, we know already such a 
state, an exact eigenfunction of zero energy, i. e. put k = 1 in (|T6|). Expanding 
the exponential in (|15|) up to first order in Zi (and substituting £* — > — 2/£) 
one finds a function of the form 

N N 

H( Zi - = J2(-0 N ~ 1 *t' MN{q)+l) (is) 

i=i q i=o q 



i. e. , the quasihole wave function earlier proposed by Laughlin [T(| as a good 
trial wave function for the CM. The r. h. s. of (|T8| ) defines the expansion in 
components with angular momenta M = M^(q) + I. Thus, there are N + 1 
degenerate states with zero energy as on the sphere. 

The inclusion of the Laughlin wave function ( |T3| ) in (0) as the quasihole 
component / = is the consequence of a somewhat ambiguous situation. The 
quasihole has a local charge deviation from the homogeneous value j- of the 
charge distribution of the electrons on the disk. E. g. for v = 1, the one-particle 



density of a quasihole at £ is [|55] 



n ( "' 5) (z) = — (l-e"5l^l 2 ) . (19) 
2n 

Around £ there is a charge depletion of magnitude e* = ve. For a finite 
system, the charge missing at £ accumulates at the finite edge reflecting the 
enlargement of the disk area by one flux quantum. The component of (|i~8|) with 
the largest possible total angular momentum M = M^{q) + iV is identical 
to the quasihole at £ = 0. If the parameter £ is considered as a particle 
coordinate for a quasihole in the lowest Landau level, smaller angular momenta 
of the electrons correspond to larger radii for the quasiholes JTI ]. Then, the 



component of ([TJD with the smallest M = M^(q) corresponds to the quasihole 
on the border of the disk and cannot be distinguished from the Laughlin state. 
In this sense, the Laughlin state (|i~3D is a quasihole too and therefore, for the 
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creation of a quasihole at a location £ 7^ 0, this is also needed. 

The angular momentum component corresponding to a quasihole at £ = 

can be chosen to generate the other quasihole components with the help of 

the magnetic translation as follows. Starting with ni=i zi^± > applying Tt, 

1 

and annihilating all total angular momentum components with a filling factor 
corresponding to m max > (N — l)q + 1 one recovers the quasihole function 
(|18|). Physically, this means a shift of the quasihole at £ = against the 
homogeneous background of the Laughlin wave function followed by a cut off 
at the border of the original system. 

Now let us turn to the spectra for one quasihole. For the Coulomb case, none 
of the spectra for 7 electrons in Fig. 2 shows any evidence of a degeneracy 
around v' = 1/3 in the low energy region. This observation expresses again 
the difficulty with the interpretation of the spectra with Coulomb interaction 
for a finite number of electrons. Therefore, coming back to the SRM, we discuss 
Fig. 4d (7 electrons). Actually, the ground state is (iV + l)-fold degenerate and 
has zero energy and angular momenta reaching from M = Mn to M^+N (here 
M = 63 to 70). This is obviously in agreement with the fact that the state 
([Dp is a zero energy eigenstate, which follows because the multiplication with 
powers of Zj can not change the minimum relative angular momentum from 
3 in Laughlin's wave function. The low energy states, appearing in branches 
above the zero ground states with excitation energies of an order of magnitude 
smaller than the energy scale given in the SRM by Vf RM ~ 0.44, are connected 
with the edge excitations U56| . 

We want to determine the quasihole energy in the SRM. Since there are three 
external parameters in our model, three mechanisms exist for the creation 
of quasiparticles. Correspondingly, the three energies (the indices ± mean 
quasihole and quasielectron, respectively) are the gross quasiparticle energy e± 
(change of the electron number), the proper quasiparticle energy e± (change 
of the magnetic field) and the neutral quasiparticle energy e\ (change of the 



area) [|27| , |55| . There are relations between these energies involving the ground 
state energy per particle of the stable state. Thus, it is sufficient to know one 
of the energies. The energy taken from our spectra is the neutral quasiparticle 
energy because we change the disk area for the creation of the quasiparticle. 
In the SRM model with zero ground state energy, all these energies are equal, 
i. e. , e = i = e n {y < 1). The neutral quasihole energy for the v = 1/g-state 
is defined as the difference of the ground state energy of the one-quasihole 
spectrum (m max = q(N—l) + l) and that of the stable state {m max = q(N—l)). 
We find from the exact spectra in agreement with the analytical arguments 
given above e1(l/3) = e_(l/3) = £-(1/3) = 0. 

Two-quasihole spectra are pictured in Fig. lc for 6 electrons and in Fig. 4e 
for 7 electrons (m max = 20), where we have again a degenerate ground state. 
But in contrast to the one-quasihole spectrum, where the zero energy levels 
in one block (a line) were non-degenerate, the levels in this case are partially 
degenerate; e. g. , in Fig. 4e we find for the M reaching from 63 to 77 (M = M 7 
to M = M 7 + 2 • 7 = 77) energy degeneracies of 1, 1, 2, 2, 3, 3, 4, 4, 4, 3, 3, 
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2, 2, 1, 1. This subspace of the lowest energy eigenstates (here of zero energy) 
can be considered as being made up of bosonic quasiholes. The quasihole can 
occupy N + 1 one-particle states, thus d bosons (quasiholes) form a Hilbert 
space of dimension (^ V+1 ^ +d_1 ) = (^J. In particular, the total degeneracy 
for two quasiholes is just (^.j -2 ), i- e. , for iV = 7, there are 36 states, and the 

exact spectrum actually shows this. In general, the degeneracy gs(N, d, M) of 
a level with M = M^(q) + M in a d-particle Bose system with N + 1 one- 
particle states can be derived from the following generating function which is 



analogous to the fermionic formula (12) 



d 1 _ T (N+l)+(l-l) Nd 

II \_ x i = E 9b(N, d, M)x M . (20) 

1=1 X M=o 

Additionally, we confirm the validity of this description by inspection of the 
degeneracy of the zero energy levels of, e. g. , the three-quasihole spectrum 
in Fig. 4f (d = 3,N + 1 = 8, M = 0, ...,21). On the other hand, it is 
immediately clear that the total dimension of the Hilbert space is actually 
Z%i o9B (N,d,M) = ( N + d ). 

We summarize that the 1/3-quasihole spectra in the SRM corraborate strongly 
the quasihole picture already for finite N. Unfortunately, the interaction be- 
tween the quasiholes is zero. A bosonic description of the quasiholes appears 
quite natural considering the degeneracies in our spectra. This conclusion is in 
agreement with investigations on the sphere where the Hilbert-space dimen- 
sion can be understood in the boson-picture E3, E4 . 



7 Quasielectrons 

In the hierarchical construction, the quasiparticles enter symmetrically, but in 
their energies and spectra, they are not very symmetric |53fl . For the energies, 
this is quite inevitable, since a 1/3-quasielectron is created by decreasing the 
maximum angular momentum to m max = q(N — 1) — 1 in order to reach a 
larger v. Confining the electrons to a smaller disk gives rise to a stronger 
interaction contribution compared to the case of an enlargement of the disk. 
Thus, (1/3) must be positive. 

Since the SRM plays the role of a canonical model, particularly with respect 
to the quasiholes, we hope that it reflects the essential physics in the quasi- 
electron case too. Still, the situation for the quasielectrons is different from 
that for the quasiholes because we can not expect zero energy eigenvalues and 
we do not know an exact eigenfunction for the quasielectron up to now. 
Let us look at the one-quasielectron spectrum for the SRM for 7 electrons in 
Fig. 4b. The ground state at M-j = 63 has a non-zero energy. This is to be 
expected because it is not possible to avoid relative angular momentum one in 
any wave function with M < and m max = 3(N — 1) — 1. Further, there are 
low lying energy levels with M > (see Fig. 4b) which are, as in the case of 
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the quasihole, again attributed to the edge of the system. The quasielectrons 
are the lowest energy levels with M = — N, ... , Mjv — 1 (here from M=56 
to 62) whose energies are almost degenerate. 

While this seems to be in contradiction to the usual definition of the quasielec- 
tron as a ground-state of the spectrum, we justify our focus on these angular 
momentum blocks by considering the quasielectron trial wave functions. That 

re- 



there are many proposals for such trial wave functions \Z7} §5], |57|, |58| 
fleets the fact that one has for the quasielectron a more complicated situation 
than for the quasihole. Still, all proposals have angular momentum compo- 
nents in the range of M — — N, . . . , M^. This becomes particularly clear, 
if on considers Haldane's quasielectron on the sphere This iV + 1-fold 

degenerate wave function is a good trial wave function for the description of 
the quasielectron which is well exhibited in the spectra on the sphere 



After a projection onto the plane according to the rules given in |28| the re- 



sulting wave function has N + 1 angular momenta components reaching from 
M = — N to Mjv and this is our trial function introduced in f57f. The 



component with M = has much smaller energy than the other components 
and is therefore excluded. 

After having identified the angular momenta of the quasielectron components 
we want to evaluate the energies. The exact energies for up to 8 electrons have 
been compared with those of Laughlin's, Jain's, MacDonald's and our quasi- 



electron trial wave functions in our previous work [3~7|. We extend here this 
comparison to 9 electrons. The exact energies for N = 9 are pictured in Fig. 
5a together with the trial function data from Table lb in P7|; the data from 



Table lc in [37| for N = 10 (without results from exact diagonalizations) are 
shown in Fig. 5b. For both particle numbers, the pictures are very similar with 
respect to the M-dependence of all trial wavef unctions; the exact eigenener- 
gies are almost independent of M. Jain's quasielectron wave function has, in 
comparison with the other trial functions, the best energy expectation values, 
but for higher M, the deviations from the exact values are still quite large. 
The exact data in the one-quasielectron spectrum exhibit already for a finite 
electron number almost degeneracy. This degeneracy is not so well expressed 
in Fig. 4b as in Fig. 5a, because of the different scale; the relative deviation 
from the average value is for N = 9 only 5.6 %. 

Next we try to find the neutral quasielectron energy £"(1/3) from the exact 
numerical data for up N = 9 electrons in the SRM. Unfortunately, an extrap- 
olation with respect to N for a single component of the quasielectron, e. g. for 
the one with the smallest angular momentum M = — N, does not work 
because the energy fluctuates with N. Therefore, we calculated the arithmetic 
mean of the quasielectron energies of all N components forming the new quan- 
tity < £"(1/3) >. Then, these data are extrapolated to the thermodynamic 
limit by minimizing the weighted quadratic deviation from an extrapolation 
curve with a finite size correction proportional to a negative rational power of 
iV to be determined. This way, we find in the thermodynamic limit 

< £+(1/3) >= 0.1865 + 1.8790 jy -2 ' 6320 (21) 
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This curve is shown in Fig. 6a. The finite size corrections vanish surprisingly 



fast. Arguments given by Morf et al. p7|, [59[| for a iV 2 correction are not 



applicable here because there is no background in the SRM. The result ( |2lD 
can be compared with the numerical data for the discontinuity in the chemical 
potential in the case of short-range interacting electrons on the surface of the 
sphere. We extrapolate from Table 2 in [^6| the value for the jump of the 
chemical potential and find that the neutral quasielectron energy at v = 1/3 
is 0.1905 showing quite good agreement between our extrapolation and the 
result of the other geometry. 

We now extrapolate the energy expectation values of the trial wave functions, 
see Fig. 6b, again using the average of the available components at a given N. 
Here, the leading finite size correction is unknown too. The extrapolation of 
Jain's quasielectron is closest to the exact data. On the other hand, the finite 
size correction is qualitatively different from that of the exact data, see Fig. 
6a, particularly the approximate 1 / v^/V-behavior of Jain's trial wave function 
is not found in the exact data. Still Jain's proposal comes closest to the 
extrapolated value of 0.1865. The proposal of jfjTj, on the other hand, yields a 
higher limit, but the behavior of its finite size correction is much more similar 
to the exact result. 

One of the most interesting questions for a test of the relevance of the notion 
of a quasielectron is the understanding of the many-quasiparticle behavior. 
The simplest case is the two-quasielectron spectrum shown for 7 electrons 
for the SRM (Fig. 4a). We realize that different from the case of the two- 
quasihole spectrum there is a complicated interaction superposed by finite size 
effects. It seems difficult to identify the two-particle states arising from the 
subspace generated by the quasielectrons. One of the reasons should be a 
strong interaction between the quasielectrons leading to a mixing with other 
energy levels of the many-electron system. Thus, a detailed investigation of 
the quasielectron interaction for two quasielectrons and further implications 
for more quasielectrons in the framework of the hierarchical scheme remain a 
matter of future research. Now, there is an intensive discussion on this point 



17, 54, 50 1 and a satisfying explanation of the success of Jain's wave function 



and its connection with the ideas of the hierarchical scheme are still missing. 



We do not believe that a proposal reconciling the two pictures [fKj is justified 
up to now. 



8 Summary 

The numerical diagonalization of the FQHE Hamiltonian for up to nine elec- 
trons on the disk shows that whether a unique identification of a particular, 
stable state at filling factor 1/q is possible in the spectra depends on the elec- 
tronic interaction. This is in contrast to the spherical geometry. Whereas 
for the Coulomb interaction, the ground states of spectra around 1/3 can not 
be related to the Laughlin wave function for small N, the latter turns out to 
be the exact ground state solution of the model with short-range interaction. 
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The picture that the spectra for Coulomb and for the short-range interaction 
are very similar and that the states of both interactions are related to each 
other by only a small perturbation could not be quantitatively confirmed for 
the present geometry of a finite disk. 

For the identification of the stable state and the quasiholes, the SRM serves as 
a canonical model. Then, N+l zero-energy levels can be exactly identified in 
the one-quasihole spectra as quasihole components of the Laughlin type. The 
degeneracy of these levels in the spectra for two and more quasiholes justifies 
to treat the quasiholes as bosons. 

In the quasielectron case, N components of the quasielectron are identified 
which show almost degeneracy with small deviations from the average value in 
our finite size calculations. The expectations values of trial wave functions of 
various authors are compared with the exact data. They deviate with increas- 
ing angular momentum more and more from the exact values leaving space 
for improvement. An extrapolation of the average quasielectron energy to the 
thermodynamic limit is in agreement with the results of other geometries. 
It should be interesting to use the quasiparticle identification on the disk in 
order to extend the present treatment to the quasielectron-quasielectron in- 
teraction. Creation energies and interaction parameters could then be used 
as input parameters in an attempt to formulate an effective theory for the 
quasielectron. 
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Table 1 



N 


iV$ = m max + 1 


Mg 


Remarks 


5 


13 
14 
15 
16 
17 
18 
19 


25 
25 
30 
30 
30 
30 
35 


v = 1/3 

v — 1/3 + one quasihole 

v — 1/3 + two quasiholes, u' — 1/3 




13 


30 


Laughlin state 


6 


16 
17 
18 
19 
20 
21 
22 


35 
39 
45 
45 
45 
45 
51 


v = 1/3 

v — 1/3 + one quasihole 

v — 1/3 + two quasiholes, v' — 1/3 




16 


45 


Laughlin state 


7 


19 
20 
21 
22 
23 
24 
25 


51 
57 
57 
63 
63 
63 
69 


v = 1/3 

v — 1/3 + one quasihole 

v — 1/3 + two quasiholes, v' — 1/3 




19 


63 


Laughlin state 



Search for the stable state" in the Coulomb model (CM) 
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Figure Captions 



Fig. 1: Complete many-particle energy spectrum for 6 electrons and various 
interactions. The levels are ordered with respect to the total angular momen- 
tum M. 

a: Coulomb interaction with background, degeneracy 18 [ra max = 17), i. e. 
v' = 1/3. 

b: as a, but without background. 

c: as a, but for the short-range interaction (|8p without background, 
d: as c, but with background. 

e: short-range interaction without background, m max = 15, i. e. v — 1/3. 

Fig. 2: Series of spectra showing for decreasing filling factor the low lying 
states for 7 electrons with Coulomb interaction; m max = 16, . . . , 21 (Fig. 2a - 

f)- 

Fig. 3: Energy per particle for the model with Coulomb interaction and back- 
ground at v' = 1/3 in dependence on the number of electrons N up to iV = 9. 
The line extrapolates to the thermodynamic limit. 

Fig. 4: Series of spectra showing for decreasing filling factor the low lying 
states for 7 electrons with the SRM interaction (|j); m max = 16, . . . , 21 (Fig. 
4a - f). Fig. 4c shows the stable state, Fig. 4d and 4b the one-quasihole state 
and the one-quasielectron state, respectively. 



Fig. 5: Numerically calculated expectation values for Laughlin's, our ]57] , and 
Jain's quasielectron trial wave functions (from top to bottom) arranged with 
increasing M; 

a: iV = 9 electrons, the lowest levels (dottet line) are the exact energies, 
b: N = 10 electrons. 

Fig. 6: Averaged quasielectron energy in the SRM @ extrapolated to the 

thermodynamic limit. 

a: exact energies for up to N = 9. 

b: Jain's (x), Laughlin's (+) and our (*) ]57| quasielectron trial wave functions 
for up to iV = 10. 
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